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Abstract — The voltage-carrying state of superconducting 
NbTiN nanowires, used for single-photon detectors, is analyzed. 
Upon lowering the current, the wire returns to the supercon- 
ducting state in a steplike pattern, which differs from sample 
to sample. Elimination of geometrical inhomogeneities, such as 
sharp corners, does not remove these steplike features. They 
appear to be intrinsic to the material. Since the material is 
strongly disordered, electronic inhomogeneities are considered as 
a possible cause. A thermal model, taking into account random 
variations of the electronic properties along the wire, is used as 
an interpretative framework. 



I. Introduction 

Superconducting nanowire single-photon detectors 
(SNSPD's) are promising devices for the detection of 
single photons because of their fast response time, broadband 
sensitivity, and low dark-count rate JT] . A variety of materials 
are currently being used, with the common denominator 
that their resistivity in the normal state is unusually large, 
consistent with an electronic mean free path in the order of 
the interatomic distance. These materials are chosen for their 
high critical temperature and fast electron-phonon scattering. 
In view of the increased interest, the need for reproducible 
device-fabrication has become an issue. It has been found 
that often variations of detection efficiency occur from device 
to device 12 and as a function of position on a single device 
|3l , which are not easily understood. Part of this lack of 
control may signal the need for improved materials control 
on an atomic level. However, since SNSPD's are made from 
highly resistive superconductors, an intrinsic cause might be 
present as well. 

Due to the short elastic-scattering length, localization is 
competing with superconductivity, and a tendency to an insu- 
lating state is accompanied by a tendency to become supercon- 
ducting. It is predicted that intrinsic electronic inhomogeneities 
are formed irrespective of the specific atomic inhomogeneity 
10) . Recent experimental work on TiN has clearly demon- 
strated the occurrence of such electronic inhomogeneities by 
measuring the local values of the superconducting energy- 
gap by scanning tunneling microscopy ]5). Stimulated by 
these observations, we have recently studied the microwave 
electrodynamics of such films (6), and analyzed the data in 
the context of a recent theory proposed by Feigel'man and 
Skvortsov 0. Unfortunately, for SNSPD's a limited number 
of parameters is available to characterize them: the critical 
temperature, the critical current, the resistivity and the resistive 
transition. It is, however, reasonable to assume that variations 
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Fig. 1 . Examples of the steplike decrease in voltage observed for decreasing 
current-bias for various wire geometries. The inset shows the temperature 
dependence of the resistance of one of the 20 /^m-long, 1 /^m-wide, and 
8 nm-thick wires, including the gradually widening contact pads (§). 



of the superconducting gap, as observed by Sacepe et al. play 
a role in the observed variation of detection efficiency from 
device to device and from position to position on the device. 
In this manuscript, we focus on one of the few available 
extra sources of information; the return of the device to the 
superconducting state, the so-called retrapping current (Fig. [T}, 
which is known to vary from sample to sample and shows a 
clear steplike structure. 

II. Experiments 

The measurements described in this paper, were performed 
on the same NbTiN wires described in [8||, where the reduction 
of the critical current due to current crowding was studied. 
The NbTiN films have a thickness d of ~8 nm, a critical 
temperature T c of 9.5 K, and a resistivity p of 170 /ii7cm. The 
wires have a width w of 1 /im. We have chosen this width, 
though it is wider than typical SNSPD devices, since wider 
wires are less sensitive to fabrication defects and show stronger 
current-crowding effects. Similar results were obtained for 
samples as narrow as 50 nm. The length L is 20 /im. Four 
different wire geometries were studied: straight wires, corners 
with optimally designed inner curves, corners with sharp inner 
curves, and a straight line that included a stub extending 1 /im 
to one side. 

For each of the devices the current- voltage characteristics 
were measured in a 4.2 K liquid helium dipstick. The samples 
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were mounted on a ceramic PCB that was in direct contact 
with a copper sample stage with thermally anchored wires, 
to ensure a temperature equal to the sample-stage. The tem- 
perature of the sample stage itself was controlled by a heater 
element and measured by a temperature sensor embedded in 
the sample stage. The measurements reported here are all 
current-biased in a two-point measurement configuration. The 
wires were driven into the normal state by applying a current 
above the critical current of the wire. Then we monitored for 
each sample for decreasing current the return from the normal 
state into the superconducting state. 

III. Steplike transitions 

Fig. Q] shows the retrapping characteristics for four ge- 
ometries. This set of measurements was taken at 4.2 K. In 
all samples a stepwise transition from the normal state to 
the superconducting state is observed. The intermediate states 
between the normal and superconducting state of the wire are 
stable for long periods of time (measured up to ~10 minutes) if 
the bias current is kept stable at this point. The specific pattern 
of a wire is reproducible over multiple measurements. We 
observe similar stepwise retrapping characteristics in 100 nm 
wide NbTiN and TiN nanowires. We find that the presence and 
stability of these steps does not depend on the geometrical 
shape of the wire. All geometries show stepwise retrapping 
patterns. This shows that, for these geometries, the presence of 
a geometrical 'constriction' does not dominate the retrapping 
characteristics. This is clearly different from the supercon- 
ducting to normal-state transition, where the critical current is 
significantly suppressed in geometries with sharp features (8|. 
The available evidence suggests that these disordered super- 
conducting wires have a tendency to become superconducting 
by forming one or more superconducting domains, which 
expand in a steplike discontinuous way. 

We have recently performed a detailed study of a model 
system consisting of a superconducting wire between two 
normal reservoirs. The analysis was executed taking the full 
non-equilibrium processes into account J9). In the strongly 
disordered wires of NbTiN (and TiN) the origin of the steps 
needs to be determined prior to a more detailed analysis. It 
has been theoretically suggested that for a superconductor 
with homogeneous disorder the superconducting state tends 
to develop mesoscopic fluctuations in the value of the energy 
gap and a flattening of the peak in the density of states 
J3). Scanning tunneling microscopy measurements on thin 
films of TiN, with a resistivity even higher than our films: 
p > 1000/iSlcm, show such a variation of the gap energy 
over the film Q. The maximum variation of the gap energy 
is about 20% of the total gap energy, and the typical size of 
the regions is about 50 nm, much larger than the grain size of 
the films. Variations in the geometry of the nanowire, such as 
thickness variations due to steps in the underlying substrate, 
might lead to similar behavior as modeled here [TTOl . We use 
SEM-inspection to guaranty the homogeneity of the width. We 
observe similar behavior in ALD-deposited TiN nanowires on 
an amorphous silicon oxide substrate. The persistent presence 
of these steplike features for wires that appear 'perfect', leads 
us to conjecture that electronic inhomogeneities play a role. 
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Fig. 2. Schematic representation of the model, (a) To model inhomogeneities, 
one can let the critical temperature, heat transfer coefficients, or resistivity 
vary as a function of position in the wire, (b) In each slice of the wire, the 
electron and phonon system are assumed to have a well-defined temperature. 
Energy enters the electron system due to Ohmic dissipation in the normal 
regions. Heat diffusion along the wire occurs both in the electron and in the 
phonon system. Heat is exchanged between the electrons and the phonons, 
and between the phonons and the substrate. 



IV. Thermal model with randomly fluctuating 

PARAMETERS 

We assume that the transition occurs by the formation of 
expanding superconducting domains, with the observed resis- 
tance given by the backscattering resistance in the remaining 
normal part of the wire. In principle, for a driven system 
consisting of superconducting and normal domains one would 
have to include current conversion processes at the interfaces. 
Since we are interested in the possible influence of random 
electronic inhomogeneities, we use the much more tractable 
thermal model of the superconducting wire introduced by 
Skocpol, Beasley, and Tinkham ifTTI . 

In Fig. |2] we illustrate our model. The bulk contacts are 
assumed to be in the superconducting state and both the 
electron- and phonon-temperature are taken to be at the bath 
temperature of 4.2 K at the contacts (at x = and x = L in 
Fig. |2 a)). The heat flow model in each segment of length dx is 
shown in Fig.|2jb). We assume that in each slice of the wire the 
electrons and phonons each have a thermal distribution, with 
temperatures T c and T p respectively. The wire is biased with 
a constant current /, and has a normal state resistance R n . We 
also assume the temperature is constant over the width of the 
wire. The electron system looses heat by diffusion along the 
wire and by electron-phonon interaction. The phonon system 
looses energy by phonon-diffusion along the wire and escape 
of phonons from the film to the substrate. 

We write a static heat balance equation for both the electron 
and phonon baths: 

-Pdc + Pdi$,e( X ) ~ p dis,e(x + dx) 

+Pph-e - Pc-ph = (1) 

-P csc + Pdiff,ph(aO - Pdiff, P h(£ + dx)+ 

Pe-ph - Pph-e = (2) 

We describe the heat transport by electron and phonon 
diffusion with Fourier's law: 



Pdiff.o = -\ c —wd, 
dx 

dTp 

"diff.ph = —Xph—j^-wd, 



(3) 
(4) 
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with wd the wire cross section, and A c and A p h the electron 
and phonon heat transfer coefficients: 

2 

aT e , (5) 



x 



ph 



(6) 



with a the normal state conductivity, c p h the phonon heat 
capacity, and _D p h the phonon diffusion coefficient. Assuming 
two-dimensional phonons, the phonon heat capacity is given 

by 

Cph « 43.3A r modcs fcB (J!j~j ' ( 7 ) 

with iV m odcs the number of phonon modes per unit volume, 
and Od the Debye temperature. The phonon diffusion coeffi- 
cient is given by 



D 



1 



P h 



(8) 



with u the speed of sound in the material, and d the film 
thickness. The electron-phonon coupling is parameterized by 
a time constant T _ p h. The energy flows are given by 



Pc 



- P h 



P 



P h- 



7~e— ph 
Cph 

7ph— c 



w 



d ■ dx, 



d ■ dx, 



(9) 
(10) 



with c c = \{^f-) 2 jy-T e the electronic heat capacity, esti- 
mated with the Fermi free electron gas model, with D c the 
electronic diffusion constant. Phonon escape to the substrate, 
which is assumed to be at the bath temperature, is also 
parameterized by an escape time t csc : 

Cph 



P. 



(T p — Tb) w ■ d ■ dx. 



(11) 



In this straightforward model, the crucial part is the electron 
temperature and electronic heat flow. Assuming that electronic 
inhomogeneities play a role, it may lead to random variations 
in T c . However, we find that the resistive transition R(T) 
curves of the wires are very sharp (inset Fig. [Q. An alternative 
source of fluctuations might be contained in the electronic 
heat conductivity, which would then also be present in the 
resistivity. In principle, the superconducting state is insensitive 
to elastic scattering. Variations in elastic scattering can easily 
be tolerated, while maintaining a uniform T c ifTZl . However, in 
these strongly disordered films the superconducting properties 
do depend on disorder and a variation of T c with position is 
possible J6). 

For simplicity, we have chosen to model the intrinsic 
electronic inhomogeneities by letting the critical temperature 
vary randomly along the length of the wire. At points at which 
the electron temperature exceeds the local critical temperature, 
the wire is considered to be in the normal state and Ohmic dis- 
sipation is present. At points at which the electron temperature 
is lower than the critical temperature, the wire is considered to 
be in the superconducting state and no dissipation is present: 



Pdc = 



I 2 R n f, if T e (x) > T c (x) 
0, if T c (x) < T c (x) 



(12) 
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Fig. 3. Electron temperature profiles along a wire for decreasing current- 
bias from 122 (red) to 61 fiA (blue). The black line represents the randomly 
varying critical temperature along the wire. The stepwise transition from 
normal state to superconducting state can be seen from the temperature profiles 
in which only a part of the wire has an electron-temperature above the critical 
temperature. 



Several parameters play a role but many are known by 
design (the geometric parameters w, d, and L), or can be mea- 
sured independently (D e , R n , a). The following parameters 
were estimated from literature values for thin NbN films and 
applied to our NbTiN films: 

• Tc-ph- We use an empirical relation T c _ p h « 500^ 16 
ps, that was reported for NbN thin films IfPJl . We use the 
value at T c and ignore the slight temperature dependence. 

• r osc . The value reported for 3.5 nm NbN on sapphire is 

t csc = 38 ps m. 

• N mo d cs - This value is taken equal to the atomic density, 



for NbN this is 



10 



28 



na. 



• 0d- The two-dimensional Debye temperature is Q^f w 
0.919^. Ql The value of the latter is 9^° = 250 K 

ma. 

. u = 2.3 km/s for NbN ED- 

V. Simulated characteristics 

In Fig. [3] we show two sets of electron temperature profiles 
for a 20 /im long, 1 /im wide, 8 nm thick, NbTiN wire, 
simulated using the model described in section [IV] The black 
line represents the critical temperature, which varies randomly 
around the experimentally determined value T c ,avg = 9.5 K. 
The used parameters are AT C = 0.5 K, where T c {x) = 
7c,avg + x- AT C , with x a random number between — ^ and i. 
The length of the T c regions of the random pattern is 500 nm, 
a value arbitrarily chosen but suggested by the experiments by 
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Sacepe et al. (5). In the two panels of Fig.[3]a different random 
pattern with the same parameters is chosen, clearly leading to 
a different approach to the superconducting state. The bias 
current is decreased from 122 /xA, at which the entire wire is 
normal. As the current is lowered, the electron temperature in 
the wire decreases smoothly, until a stepwise jump occurs in 
which the superconducting regions at the boundary increase 
to a new stable position. This continues until finally the entire 
wire is in the superconducting state. 

The phonon-temperature follows the same profile as the 
electron-temperature. In the normal region it lies between 5.9 
and 6.3 K for the currents shown. Because of the low electronic 
diffusion constant for our films, D « 1 cm 2 /s, and strong 
electron-phonon coupling, the energy flow is dominated by 
energy transfer to the substrate via the phonon system. In 
other words, the size of the normal regions is always much 
larger than the thermal healing length in the wire. Therefore, 
under current-bias, the electron- and phonon-temperature s are 
constant over the normal region and their value is independent 
of the size of the normal region. 

Fig. |4] shows the current-voltage characteristics which re- 
sults from the electron-temperature profiles in Fig. |3] Elements 
of the wire where the electron temperature is above the 
critical temperature are taken to be in the normal state, with 
a resistance, whereas elements of the wire where the electron 
temperature is below the critical temperature are taken to be 
in the superconducting state. Jumps in the transition from the 
normal state to the superconducting state can be seen around 
95 /iA. The precise value of the retrapping current differs from 
the experimentally observed value of ~ 70 /iA. However, we 
are not aiming for a quantitatively accurate prediction of the 
retrapping current, since a number of thermal parameters are 
not well-known. We are aiming for possible contributions to 
the observed steplike features. The exact shape depends on the 
particular random inhomogeneity pattern of the sample. 

Steps in the current-voltage characteristics are also reported 
in ifTTl . where they are ascribed to phase slip centers. These 
current-voltage characteristics are however qualitatively differ- 
ent from the ones we observe. In IfTTl , the stepwise transition 
occurs mainly in the increasing current branch, in which we 
observe an immediate transition from the superconducting- to 
the normal-state. 

VI. Conclusions 

We observe a persistent stepwise transition from the normal 
to the superconducting state under decreasing current-bias in 
NbTiN wires, for all wire geometries. It is argued that the 
stepwise pattern may be due to disorder fluctuations or elec- 
tronic inhomogeneities intrinsic to disordered superconductors, 
modeled by local variations in the critical temperature. 

Acknowledgment 

The authors would like to thank Karl Berggren for stimulat- 
ing discussions. We acknowledge support from the Foundation 
for Fundamental Research on Matter (FOM), Microkelvin 
(No. 228464, Capacities Specific Programme), and the Teradec 
program of NWO-RFBR. 



1 1 1 
-•-Fig. 3(a) 


i 1 i 1 


-•-Fig. 3(b) 






|t» 




1/ 


, i 


T 

i.i. 



85 90 95 100 105 

I(mA) 

Fig. 4. Current-voltage characteristics deduced from the simulated electron 
temperature profiles in Fig. [3] Steps in the transition from the normal to the 
superconducting state are seen for decreasing current-bias. 
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